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ABSTRACT 



Aims. We study the influence of rotation and disc lifetime on lithium depletion of pre-main sequence (PMS) solar-type stars. 
Methods. The impact of rotational mixing and of the hydrostatic effects of rotation on lithium abundances are investigated by com- 
puting non-rotating and rotating PMS models that include a comprehensive treatment of shellular rotation. The influence of the disc 
lifetime is then studied by comparing the lithium content of PMS rotating models experiencing different durations of the disc-locking 
phase between 3 and 9 Myr. 

Results. The surface lithium abundance at the end of the PMS is decreased when rotational effects are included. During the beginning 
of the lithium depletion phase, only hydrostatic effects of rotation are at work. This results in a decrease in the lithium depletion 
rate for rotating models compared to non-rotating ones. When the convective envelope recedes from the stellar centre, rotational 
mixing begins to play an important role due to differential rotation near the bottom of the convective envelope. This mixing results 
in a decrease in the surface lithium abundance with a limited contribution from hydrostatic effects of rotation, which favours lithium 
depletion during the second part of the PMS evolution. The impact of rotation on PMS lithium depletion is also found to be sensitive 
to the duration of the disc-locking phase. When the disc lifetime increases, the PMS lithium abundance of a solar-type star decreases 
owing to the higher efficiency of rotational mixing in the radiative zone. A relationship between the surface rotation and lithium abun- 
dance at the end of the PMS is then obtained: slow rotators on the zero-age main sequence are predicted to be more lithium-depleted 
than fast rotators due to the increase in the disc lifetime. 
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1. Introduction 



Rotation is one of the ke y processes to influence all outputs 
of stellar models (see e.g. 'Maeder"2009'). While rotational ef- 
fects have an especially strong impact on the physics and evolu- 
tion of massive stars, they can also change the global properties 
and evolution of solar-type stars by counteracting the efFets of 
atomic diffusion an d bringing fresh hydrogen fuel to the cen- 
tral stellar core (e.g. Eggenberger et al. 2010b). Observations of 
light element abundances bring valuable constrain ts to progress 
in our modelling of these tra nsport processes (e.g.lPinsonneaulti 
I2OIOI ). In a preceding work (Eggenberge r et all 2010aV we in- 
vestigated the influence of shellular rotation on the lithium abun- 
dances of main-sequence solar-type stars in the context of a pos- 
sible link between lithium depletion and the r otatio nal history 
of exoplanet-host stars proposed by iBouvie idlooi. We then 
found that the larger efficiency of rotational mixing predicted 
in exoplanet-host stars can result in lower lithium abundances 
for these stars during the main sequence. After this study of the 
impact of shellular rotation on the lithium abundances of main- 
sequence stars, we are now interested in investigating the effects 
of rotation on the lithium content of pre-main sequence (PMS) 
solar-type stars. 

The impact of rotation on PMS lithium depletion o f solar- 
type s tars has alr eady been discussed (e.g. Pln sonnea ult et al.l 
119901: lMartin& Cl aret 1996t iMendes et al.lTl999l) . These previ- 
ous studies focused mainly on the hydrostatic effects of rota- 



tioifl and did not consider the influence of a possible coupling 
between the star and its disc. The star-disc interaction plays a 
key role in determining the rotation of the star on the zero-age 
main sequence (ZAMS). By assuming that the presence of the 
disc prevents the star from spinning up so that its surface an- 
gular velocity remains constant, one finds that the influence of 
this coupling simply relies on the disc hfetime. A longer disc 
lifetime then enables the star to lose a larger amount of angular 
momentum during the PMS, resulting in lower rotation rates on 
the ZAMS. In addition to its effect on the surface rotational ve- 
locity during the PMS, the duration of the disc-locking phase can 
change the internal rotation profiles, with a possible impact on 
the lithium content. In this paper, we thus investigate the influ- 
ence of rotation and disc lifetime on the PMS lithium depletion 
of solar-type stars. Both effects of rotational mixing and of the 
centrifugal force are studied by computing models that include 
a comprehensive treatment of shellular rotation. The effects of 
the disc lifetime on the PMS lithium depletion are then investi- 
gated by comparing the PMS evolution of models with different 
durations of the disc-locking phase. 

The modelling of rotation is first briefly described in Sect.|2] 
The effects of rotation on PMS lithium depletion are discussed 
in Sect. [3] while the impact of the disc lifetime is investigated in 
Sect.|4] The conclusion is given in Sect. |5] 



' By hydrostatic effects of rotation, we mean here the effects related 
to the change in the effective gravity due to rotation (the effective grav- 
ity being defined as the classical gravity diminished by the centrifugal 
acceleration). 
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2. Modelling rotation 

The stellar evolution code used for these computations is the 
Geneva code, which includes a detailed treatment of shellular 
rotation ( Eggenberger et al. 2008). In the context of shellular 
rotation, the transport of angular momentum in radiative zones 
obeys an advection-diffusion equation written in Lagrangian co- 
ordinates (IZahnlll992HMaeder & Zahnl[T998l) : 



d/,\ I d I , ^ Id ,dQ.\ 



(1) 



where r is a characteristic radius, p the mean density on an iso- 
bar, and Q(r) the mean angular velocity at level r. The vertical 
component u{r, 9) of the velocity of the meridional circulation at 
a distance r to the centre and at a colatitude 6 can be written 



u(r, 0) = f/(r)P2(cos 0) . 



(2) 



Only the radial te rm Ujr) appears in Eq. ([T). Its expression is 
given by (Maeder ^Zahni[T 998) 



U{r) 



{En + Eu) 



(3) 



where P is the pressure, Cp the specific heat. En and E^^ are 
terms depending on the Q- and /i-distributions respectively, up to 
the third- order derivatives and o n various thermodynamic quan- 
tities (see Maeder & Zahn 1998, for more details). 

Meridional circulation and shear mixing are considered as 
the main mixing mechanisms in radiative zones. In convective 
zones, solid-body rotation is assumed as indicated by the solar 
case. The first term on the right-hand side of Eq. ([TJ describes 
the advection of angular momentum by meridional circulation, 
while the second term accounts for the transport of angular mo- 
mentum by shear instability with D - Dshear- The expression of 
this diffusion coefficient is given by 



shear 



[lVi + i) + (v.d-v,,d)] 



9n „dlnQ 



32 dlnr 



(4) 



with K the thermal diffusivity, /energ the fraction of the excess 
energy in the shear that contributes to mixing (here taken equal 
to 1) (iMaeder & Mevnet,2001,) . and the diffusion coefficient 
associated to horizontal tu rbulence. The usual expression for this 
coefficient is, according to lZahnl(ll992h . 



Dh = -r\2V(r) - at/(r)| 



(5) 



where U(r) is the vertical component of the meridional circu- 
lation velocity, V{r) the horizontal component, Ch a constant of 
order unity (here taken equal to 1) and a - The full 

solution of Eq. ([T]i taking U{r) and D into account gives the non- 
stationary solution of the problem. We recall here that merid- 
ional circulation is treated as a truly advective process in the 
Geneva evolution code. 

The vertical transport of chemicals through the combined ac- 
tion of vertical advection and strong h orizontal diffusion can b e 
described as a pure diffusive process dChabove r&Zahn! ll992h . 
The advective transport is thus replaced by a diffusive term, with 
an effective diffusion coefficient 



£>eff = 



\rU(r)? 
30Z)h 



(6) 



The vertical transport of chemical elements then obeys a diffu- 
sion equation that, in addition to this macroscopic transport, also 
accounts for (vertical) turbulent transport; 



df pr^ dr 



2 dXi 

r p(Deff -l-I>shear)-;!— 

ar 



(7) 



where X, is the mass fraction of element and X, represents the 
variations in chemical composition due to nuclear reactions. 



3. Effects of rotation 

To study the effects of rotation on the lithium abundances 
of PMS solar-type stars, we compute the PMS evolution of 
1 Mq models with a solar chemical composition as given by 
Grevesse & Noels (1993) and a solar-calibrated value for the 
mixing-length parameter. One model is first computed without 
rotation. The corresponding rotating model shares exactly the 
same initial parameters except for including a comprehensive 
treatment of shellular rotation. This model is computed with an 
initial angular velocity of 20 and a disc lifetime of 6 Myr 
During the disc-locking phase, the surface angular velocity of 
the star is simply assumed to remain constant. 

The evolutionary tracks in the HR diagram are shown for 
both models in Fig. [T] Including rotation shifts the track to 
slightly lower effective temperatures and luminosities as previ- 
ously found by many studies (see e.g. Pinsonneault et al.ll 19891 : 
lMartin&Cl'^ll996l: iMendes et al.l ll999'). Figure l2l compares 
the evolution of the surface lithium abundance during the PMS 
for the models computed with and without rotation starting with 
a lithium abundance A(Li) = \og[N{U)IN{H)] + 12 = 3.26 
dAsplund et alj|2009l) . At the end of the PMS, the rotating model 
exhibits a lower surface lithium abundance than the non-rotating 
one (difference of about 0. 1 dex for a total lithium depletion of 
about 1 dex). We thus see that the inclusion of rotation results in 
a global increase of the lithium depletion during the PMS. 

To investigate whether these differences in the lithium con- 
tent between rotating and non-rotating models are mainly due 
to rotational mixing or to hydrostatic effects of rotation, another 
rotating model of 1 Mq is computed by only accounting for the 
effects of rotational mixing. The initial rotational velocity and 
disc lifetime of this model are thus identical as for the rotating 
model including the full treatment of shellular rotation, but the 
effects of the centrifugal force on the stellar structure are not 
taken into account. The evolution of this model in the HR di- 
agram is also shown in Fig. [T] but the corresponding track is 
superimposed on the one of the non-rotating model. The change 
in the evolutionary tracks in the HR diagram induced by rota- 
tion for PMS solar-type stars is thus mainly due to the effects of 
the centrifugal force and not to a change in the global properties 
induced by rotational mixing. 

Concerning the evolution of the lithium content during the 
PMS, the situation is quite different from the changes observed 
in the HR diagram. Figure|2]indeed shows that the lithium deple- 
tion is greater for the rotating model including only the effects 
of rotational mixing than for the non-rotating model. While the 
effects of rotation on the PMS evolutionary tracks in the HR di- 
agram are mainly due to the centrifugal force, rotational mixing 
is found to play an important role in the increase of lithium de- 
pletion for rotating models during the PMS. As seen in Fig.|2] 
both rotating models exhibit lower lithium surface abundances 
than the non-rotating model on the ZAMS. We also note that the 
lithium depletion occurs earlier for the rotating model, when ac- 
counting only for the effects of rotational mixing, than for the 
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without rotation 
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log(T,„ [K]) 

Fig. 1. Effects of rotation on the PMS evolutionary tracks in the 
HR diagram for 1 M© models. The dotted and continuous lines 
correspond to a non-rotating model and a rotating model com- 
puted with a full treatment of shellular rotation, respectively. The 
track corresponding to a rotating model including only the ef- 
fects of rotational mixing is shown by a dashed line and is super- 
imposed on the track of the non-rotating model. The tracks end 
when the ZAMS is reached; the end of the disc-locking phase is 
indicated by a triangle. 



model that includes the full treatment of rotational effects. This 
is related to the hydrostatics effects of rotation. Rotational mix- 
ing has a negligible impact on the global stellar structure dur- 
ing the PMS, so that the temperature profile, the depth of the 
convective envelope, and the evolutionary tracks in the HR di- 
agram are identical for the non-rotating model and the rotating 
model computed without the hydrostatics effects of rotation. For 
the model including a complete treatment of rotation, the tem- 
perature at the base of the convective envelope is lower during 
the initial phase when the star is still fully convective. This is 
caused by the effective gravity of the star being reduced by the 
centrifugal acceleration term, and this results in a lower lithium 
depletion rate. Including the centrifugal force therefore leads to 
a higher lithium content during the beginning of the lithium de- 
pletion phase on the PMS. 

After this initial phase, a radiative core begins to develop at 
the centre of the star The temperature at the base of the con- 
vective envelope is slightly higher for the model that includes a 
complete treatment of rotation than for the non-rotating one dur- 
ing the end of the PMS evolution. After decreasing the lithium 
depletion rate during the beginning of the PMS, the hydrostatic 
effects of rotation are thus found to favour lithium depletion dur- 
in g the second part of th e PMS, in good agreement with the study 
bv lMendes et alJ(ll999l) . More importantly, rotational mixing be- 
gins to play a dominant role during this phase. As can be seen 
in Fig. |2]by comparing the non-rotating model to the model that 
includes rotational mixing only, the impact of rotational mixing 
on the surface lithium abundance begins to be visible after about 
3 Myr. By transporting lithium to deeper layers where it is easily 




log(age [yr]) 

Fig. 2. Top: effects of rotation on lithium depletion during the 
PMS for 1 Mq models. The non-rotating model is indicated by 
the dotted line, while the dashed and continuous lines corre- 
spond to rotating models that only include the effects of rota- 
tional mixing and a full treatment of shellular rotation, respec- 
tively. These rotating models are computed with the same ini- 
tial angular velocity and disc lifetime. The tracks end when the 
ZAMS is reached and a triangle indicates the end of the disc- 
locking phase. Bottom: location of the mass coordinate of the 
base of the convective envelope. The track corresponding to a 
rotating model that only includes the effects of rotational mixing 
is superimposed on the track of the non-rotating model. 



destroyed, such mixing generates an increase in the lithium de- 
pletion. In the present PMS models, shear mixing plays the dom- 
inant role for lithium depletion with a negligible contribution 
from meridional currents. The decrease in the surface lithium 
abundance by rotational mixing is thus directly related to the 
increase of differential rotation in the stellar interior during the 
PMS evolution. 

Figure [3] shows the internal rotation profile at different ages 
for the rotating model including a complete treatment of rota- 
tional effects. At the beginning of the PMS evolution, the star is 
fully convective and rotates as a solid body. As soon as the ra- 
diative core develops, differential rotation takes place in the ra- 
diative zone mainly as a result of the star contracting. During the 
disc lifetime, the star is assumed to be locked to its disc with a 
constant surface angular velocity. As soon as the disc disappears, 
the angular velocity at the stellar surface increases as shown in 
Fig. |3] For the present model computed with a disc lifetime of 
6 Myr, the angular velocity at the surface thus remains constant 
during the first 6 Myr Differential rotation in the radiative zone 
is enhanced by disc locking and leads to a steep rotation profile 
at the base of the convective envelope. This generates turbulence 
below the convective envelope through shear instability, and the 
resulting mixing can then transport lithium to deeper and hotter 
regions, where it is destroyed more efficiently. 

The quantitative impact of rotation on the PMS lithium con- 
tent of solar-type stars is of course sensitive to the initial angular 
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Fig. 3. Evolution of the internal rotation profile for a 1 Mq ro- 
tating model with a disc lifetime of 6 Myr The age in Myr is 
indicated for each rotation profile. 
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Fig. 4. Evolutionary tracks in the HR diagram for rotating PMS 
models of 1 Mq. The continuous, dashed, and dotted lines cor- 
respond to initial velocities of 20 Qq, 10 Qq, and SQq, respec- 
tively. The tracks end when the ZAMS is reached and the end of 
the disc-locking phase is indicated by a triangle. 



velocity adopted when the star is fully convective. To illustrate 
this point, two additional models of 1 Mq are computed with ini- 
tial angular velocities of 5 and 10 Qq. These models are also 
computed with a solar chemical composition, a solar-calibrated 
value for the mixing-length parameter and a disc lifetime of 




log(age [yr]) 

Fig. 5. Lithium surface abundance during the PMS for rotating 
1 Mq models. The continuous, dashed, and dotted lines corre- 
spond to initial velocities of 200©, 10 Q©, and 5 0©, respec- 
tively. The triangles indicate the end of the disc-locking phase. 
The tracks end when the ZAMS is reached. 

6 Myr Figure|4]shows the evolutionary tracks in the HR diagram 
of these models together with the 20 model discussed before. 
The evolution of the surface lithium abundance of these models 
is given in Fig. |5] The decrease in the lithium abundance at the 
end of the PMS is found to be more pronounced for higher ini- 
tial velocities. This mainly results from the increased efficiency 
of rotational mixing when the rotation rates increase. The hy- 
drostatic effects of rotation are also increased when the initial 
angular velocity is higher As shown in Fig.|4] the shift to lower 
temperatures and luminosities becomes more visible in the HR 
diagram when the initial rotational velocity increases. This leads 
to a slight decrease in the lithium depletion rate at the begin- 
ning of the PMS evolution when the initial velocity increases. 
As seen in Fig. |5] lithium depletion therefore occurs earlier for 
lower initial rotation velocities. 



4. Effects of disc lifetime 

In the preceding section, the effects of rotation on the PMS 
lithium depletion of solar-type stars have been discussed by 
comparing rotating models computed with the same disc lifetime 
of 6 Myr We are now interested in investigating the impact of the 
disc lifetime on the lithium depletion during the PMS. Rotating 
models of 1 Mq are computed for different disc lifetimes of 3, 
4, 5, 6, 7, 8, and 9 Myr All these models are computed with the 
same initial angular velocity of 20 Qq, a solar chemical compo- 
sition, and a solar-calibrated value for the mixing-length param- 
eter 

The evolution of the surface angular velocity during the PMS 
is shown in Fig. |6]for these models. The angular velocity at the 
surface remains constant during the disc-locking phase. As soon 
as the disc disappears, the surface angular velocity increases as a 
result of the stellar contraction until the ZAMS is reached. When 



4 



p. Eggenberger et al.: Rotation, disc lifetime and pre-main sequence lithium depletion 




6.5 



7.5 



Iog(age [yr]) 



Fig. 6. Evolution of the surface angular velocity for 1 M© models 
sharing the same initial velocity but computed with different disc 
lifetimes. Triangles indicate the end of the disc-locking phase for 
each model. The tracks end when the ZAMS is reached. 
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Fig. 8. Evolution of the surface lithium abundance for 1 Mq mod- 
els with the same initial velocity but different disc lifetimes. 
Triangles indicate the end of the disc-locking phase for each 
model. 
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Fig. 7. Evolutionary tracks in the HR diagram for 1 Mq mod- 
els with the same initial velocity but different disc lifetimes. 
Triangles indicate the end of the disc-locking phase for each 
model. 



the disc lifetime increases, the star keeps its inital surface angular 
velocity for a longer time and loses a larger amount of angular 
momentum. As shown in Fig. |6] the increased surface angular 
velocity then occurs later for longer disc lifetimes, resulting in 
lower rotational velocities on the ZAMS. 



Figure Q shows the evolutionary tracks in the HR diagram 
corresponding to these rotating models computed for various du- 
rations of the disc-locking phase. During the first 3 Myr, all mod- 
els follow the same track in the HR diagram because they share 
the same rotational history. Rotating models with shorter disc 
lifetimes exhibit higher rotational velocities as soon as the disc 
disappears, compared to rotating models that remain coupled to 
their disc for a longer time. Consequently, the hydrostatic effects 
of rotation are more pronounced for shorter disc lifetimes; this 
results in larger shifts to the red part of the HR diagram when 
the disc lifetime decreases. 

The evolution of the surface lithium abundance is shown in 
Fig. [8] While the effects of rotation in the HR diagram are found 
to be stronger for models with shorter disc lifetimes, we see ex- 
actly the opposite for the impact of disc lifetime on lithium de- 
pletion. Figure |8] indeed shows that lithium depletion during the 
PMS increases for longer durations of the disc-locking phase. 
Soon after the disappearance of the disc, the lithium depletion 
decreases compared to models with longer disc lifetimes. As 
a result, models computed with shorter disc lifetimes exhibit 
higher lithium abundances on the ZAMS than models experi- 
encing a longer disc-locking phase. 

Figure |9] shows the internal rotation profiles of these models 
with different disc lifetimes at an age of 9 Myr, which corre- 
sponds to the end of the disc-locking phase for the model with 
the longer disc lifetime. The gradient of angular velocity be- 
low the convective envelope increases when the disc lifetime 
increases. Rotational mixing is then more efficient in the inte- 
riors of models with longer disc lifetimes, which explains the in- 
creased lithium depletion shown in Fig. [8] We thus conclude that 
an increase of the disc lifetime leads to a decrease in the lithium 
content during the PMS evolution of a solar-type star owing to 
the larger amount of differential rotation generated during the 
disc-locking phase. 
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Fig. 9. Rotation profiles at an age of 9 Myr (corresponding to 
the disappearance of the disc for the model with the longer disc 
lifetime) for models computed with diff'erent disc lifetimes. 



The PMS lithium depletion predicted by non-rotating solar 
calibrated models is known to be too high to correctly reproduce 
the observations in open clusters (e.g. D'Antona& Mazzitelli 



19941 Ventura etal. 


19981 ISchlatd & Weiss 


Piau & Turck-Chiezd 


20021 iD'Antona & Montalban 



1999t 



Since including rotational effects results in lower lithium 
abundances for rotating models compared to non-rotating ones 
(with a higher impact for longer disc lifetimes), rotation does 
not solve this problem. One possibility for correctly reproducing 
the lithium abundances observed in open clusters is to reduce 
the value of the mixing-length parameter. The PMS evolution 
of two additional rotating models is then computed with a 
mixing-length parameter a - 0.6 q-q and two different disc 
lifetimes of 3 and 9 Myr Figure [TO] shows the evolution of the 
surface lithium abundance for these models during the PMS. 
As expected, the reduced value of the mixing-length parameter 
results in a lower PMS lithium depletion for both models. 
More importantly, we see that the increased duration of the 
disc-locking phase also leads to a significant decrease in the 
lithium abundance on the ZAMS, in good agreement with the 
results obtained above for solar-calibrated models. 

Two different parameters are found to determine the quanti- 
tative impact of rotational effects on PMS lithium depletion: the 
initial rotational velocity and the disc lifetime. For a given disc 
lifetime but different initial velocities, a direct correlation be- 
tween the surface velocity and the lithium abundance at the end 
of the PMS is then obtained so that fast rotators are predicted 
to be more lithium depleted than slow rotators. For a given ini- 
tial velocity but different disc lifetimes, one expects that slow 
rotators on the ZAMS are more lithium depleted than fast ro- 
tators. This is illustrated in Fig. [TT] which shows the lithium 
depletion on the ZAMS as a function of the disc lifetime for 
rotating models computed with different initial velocities. We 
see that models with longer disc lifetimes, hence lower rota- 
tion rates on the ZAMS, exhibit higher lithium depletions than 
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Fig. 10. Evolution of the surface lithium abundance for two 1 M© 
models computed with a mixing-length parameter a - 0.6 a©. 
The continuous and dotted lines correspond to disc lifetimes of 
9 and 3 Myr, respectively. Triangles indicate the end of the disc- 
locking phase for each model. 
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Fig. 11. Surface lithium depletion on the ZAMS for 1 Mq mod- 
els computed with different disc lifetimes (td) and initial veloci- 
ties (t^ini). A^Li is defined as the difference between the surface 
lithium abundance on the ZAMS of rotating models and the non- 
rotating one. 



models with shorter disc lifetimes that are fast rotators on the 
ZAMS. This agrees well with observations in the Pleiades of 
higher lithium ab undances fo r fast rotators than for slow ro- 
tators ( Soderblom et al. l [T99l . Moreover, longer disc lifetimes 
may favour the formation and migration of giant exoplanets, so 
that lower lithium abundance for exoplanet-host stars are also 
predicted in this context. This agrees well with the possible de- 
tection of differences in the lit hium content o f star s with and 
wit hout exoplanets repor ted by llsraelian et al.l (12009.) (but see 
also lBaumann et al]|2010i for a contrary view). 

We may wonder whether the results obtained here in 
the context of PMS models including shellular rotation are 
purely academic, given the need of another mechanism for 
the transport of angular momentum to account for the so- 
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lar rotation profile for ins t ance ( e.g. Mestel & WeissI 
Charbonneau & MacG regorl Il993t lEggenber ger et alJ 
Denissenkov 2010). The rotational periods of young solar-type 
stars suggest that slow rotators develop a high degree of differen- 
tial rotation between the radiativ e core and the convec tive enve- 
lope (see e.g. Bouvier 2008; De nissenkov et al.ll2010h . This in- 
dicates that the hypothesis of shellular rotation during the PMS 
evolution is perfectly valid in this case even if an additional pro- 
cess can eventually intervene later on during the main-sequence 
evolution, leading to a flat rotation profile at the solar age. 

For fast rotators, the situation is somewhat different. The ob- 
servation of rotational periods of young solar-type stars suggests 
that an efficient mechanism for the transport of angular momen- 
tum is needed for fast rotator s to correctly reproduce their spin- 
down in open clusters (e.g. iBouvie iDenissenkov et al.l 
l2010h . The hypothesis of shellular rotation during the PMS thus 
seems more questionable in the case of fast rotators than for slow 
rotators. This hypothesis is however not incompatible with the 
observation of rotational periods for fast rotating solar-type stars, 
because the strong internal coupling suggested by these observa- 
tions is not needed before the beginning of the main sequence. 
This means that a fast rotator can exhibit differential rotation 
during the PMS with an additional mechanism for the transport 
of angular momentum at work with a short coupling timescale, 
leading to an approximate flat rotation profile at the end of the 
PMS. Moreover, Fig. [8] shows that the effects of rotation and 
disc lifetime on lithium depletion occur mainly during the first 
lOMyr of the PMS evolution. The results obtained in this study 
will then remain qualitatively valid even in the case of a fast rota- 
tor that only exhibits differential rotation during the first lOMyr 
and then rotates approximately as a solid-body during the rest 
of its PMS evolution. We thus conclude that even if the physical 
mechanisms responsible for the transport of angular momentum 
in the interiors of solar-type stars are currently not fully under- 
stood, the assumption of shellular rotation during the PMS seems 
to be reasonable according to current observational constraints. 

We note that including a mechanism that enforces solid-body 
rotation during the whole PMS evolution will change the picture. 
In the case of solar-type stars, the efliciency of rotational mixing 
is indeed strongly reduced whe n sofid-body rotatio n is imposed 
(see for more details Sect. 3 of lEggenberger et"al]r2010 b). The 
effects of rotational mixing then become negligible compared to 
hydrostatic effects of rotation. As a result, the impact of disc life- 
time on lithium depletion differs from the one described above. 
For PMS models in solid-body rotation, a longer disc lifetime re- 
duces the rotational velocity of the star, which simply results in 
a decrease in the hydrostatic effects of rotation (without any in- 
crease in the efficiency of rotational mixing), hence to an overall 
decrease in the impact of rotation on lithium depletion. The re- 
lationship between slow rotation and low lithium abundance on 
the ZAMS is then found to disappear when solid-body rotation 
is enforced during the PMS. 

5. Conclusions 

The effects of rotation on the PMS surface lithium abundance 
of solar-type stars have been studied. Accounting for rotational 
effects results in an overall decrease in the lithium content at 
the end of the PMS. The effects of rotation on the PMS lithium 
content of solar-type stars are enhanced when the initial veloc- 
ity of the star increases. Moreover, the impact of rotation on 
PMS lithium depletion is found to depend on the duration of 
the disc-locking phase, during which the surface angular veloc- 
ity of the star remains constant. When the disc lifetime increases. 
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the PMS lithium abundance of a solar-type star is indeed found 
to decrease. This is due to the increase in the efficiency of ro- 
tational mixing, which is related to the larger amount of differ- 
ential rotation generated in the radiative zone during the disc- 
locking phase. This illustrates that the efficiency of shear mixing 
is more direcly related to differential rotation in stellar interiors 
than to surface rotational velocities, as also found in the case of 
stars undergoing magnetic braking (e.g. Mevnet et al. 2011). 

A link between the surface velocity and lithium abundance at 
the end of the PMS is then obtained: slow rotators on the ZAMS 
are predicted to be more lithium depleted than fast rotators due to 
the increase in the disc lifetime. This relationship is particularly 
interesting in the context of the observation of higher lithium 
abundances for fast rota tors than for slow rotators in the Pleiades 
(ISoderblo m et al."l993') and of the possible detection of differ- 
ences in the lithium content of stars with and without exop lanets 
reported by I sraelian et alj ( 12009) ; but see also IBaumann et al.l 
42010) for a contrasting view. The exact degree of coupling be- 
tween the star and its disc may of course differ from the sim- 
ple assumption of disc-locking used here, while the transport 
of angular momentum in the radiative zone of PMS solar-type 
stars may be influenced by additional physical processes such 
as magnetic fields. We are indeed still far from having a com- 
plete understanding of the various physical processes at work 
during the evolution of solar-type stars, but this study interest- 
ingly illustrates that the impact of rotational mixing depends on 
the interaction between the star and its disc. 
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